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(54) Tapped delay line based gain flattening filter 

(57) A method of modifying the spectral distribution 
of an optical signal includes the steps of splitting the sig- 
nal among a plurality of optical paths, delaying the prop- 
agation of the signal through one or more of the optical 
paths by a time duration which results in the desired 
spectral distribution when the signals at the output ends 
of each optical path are recombined, and then recom- 
bining those outputs. An apparatus for changing the 



spectral profile of an optical signal includes an input 
• waveguide and a plurality of tapped delay lines. Each of 
the tapped delay lines has a delay element configurable 
to impart a predetermined delay to the signal propagat- 
ing down the delay line. Couplers split an optical signal, 
propagating on the input waveguide, among the tapped 
delay lines. The delays imparted to the signals on each 
line are chosen so as to effect the desired change to the 
spectral profile. 
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Description 

BACKGROUND OF THE INVENTION 

5 1. Field of the Invention 

[0001] The present invention relates generally to optical filters, and more specifically to configurable filters for spectral 
shaping. 

io 2. Technical Background 

[0002] ' Actively configurable optical filters have many potential applications in optical communications, such as tem- 
poral pulse compression, dispersion slope compensation and spectral shaping. Spectral shaping is of particular interest 
in connection with optical amplifiers. The gain profile of a typical optical amplifier may be affected by numerous factors 
is such as, for example, fiber length, fiber composition, splice losses, as well as the population inversion level of the 
dopant atoms (e.g. erbium) providing the gain. Having a relatively large optical signal at the input of the amplifier tends 
to deplete the excited electron population, resulting in a different spectral gain profile than that for a small input signal. 
Average gain also varies with input signal power; 

[0003] One approach 1o dealing with this effect is to operate the amplifier at a constant inversion population, with a 
20 fixed spectral filter that provides gain flattening over the operating wavelength range. In this approach, the inversion 
level is held fixed by selectively attenuating the input signal until the desired inversion level is achieved. However, there 
may be circumstances where operating the amplifier at a constant gain level, rather than a constant inversion level, 
would be advantageous. By utilizing an actively configurable spectral filter, an amplifier could be used in a constant- 
gain mode, rather than constant inversion level, by flattening the gain profile of the amplifier for arbitrary inversion 
25 levels. One approach for realizing such a configurable spectral filter is to spatially disperse the various wavelength 
components of the input signal via a diffraction grating, then spatially modulate the dispersed components via a spatial 
light modulator or acoustically-formed phase grating. The polarization dependence and diffractive losses of the grating 
could be a drawback in some applications, however. Prisms could be used instead of a diffraction grating to eliminate 
polarization dependence and to decrease angular dispersion. Prisms require an increase in package size, however, 
30 which makes their use undesirable in many applications. 

[0004] What is needed, then, is an actively configurable spectral filter that is polarization-independent, which can fit 
in a relatively small package, and which does not suffer from diffractive losses. 

SUMMARY OF THE INVENTION 

35 

[0005] One aspect of the invention is an apparatus for changing the spectral profile of an optical signal that includes 
a plurality of optical paths, each having associated with it a delay element which imparts a predetermined delay to a 
signal propagating through the optical path. A plurality of couplers is configured to split the optical signal among the 
plurality of optical paths according to a predetermined splitting ratio. The predetermined delays and predetermined 
40 splitting ratios are chosen so as to effect the desired change to the spectral profile. 

[0006] In another aspect, the delay elements may include portions of the waveguide which have an altered index of 
refraction. The index of refraction may, in yet another aspect, be reversibly alterable. 

[0007] In another aspect, the invention includes a method for modifying the spectral distribution of an optical signal 
which includes the steps of splitting the signal among a plurality of optical paths, each having an output end, so that 

45 a portion of the signal propagates on each of the optical paths; changing the relative phases of the portions of the 
signal propagating through the optical paths; modifying the amplitudes of the portions of the signal propagating through 
the optical paths; and combining the portions of the signal at the output ends of the optical paths. The relative phase 
changes and the amplitude modifications have magnitudes that result in a desired spectral distribution when the por- 
tions of the signal at the output ends of the optical paths are combined. 

50 [0008] In another aspect, the invention includes a method for modifying the spectral distribution of an optical signal 
which includes the steps of splitting the signal among a plurality of optical paths according to predetermined splitting 
ratios, and changing the relative phases of the portions of the signal propagating through different ones of the optical 
paths; the relative phase changes and the splitting ratios having magnitudes which result in a desired spectral distri- 
bution when the portions of the signal at the output ends of the optical paths are combined; and combining the portions 

55 of the signal at the output ends of the optical paths. 

[0009] In another aspect, the invention includes an optical communication system comprising an optical amplifier 
and a spectral filter that includes an input waveguide and a plurality of tapped delay lines. Each of the tapped delay 
lines has a delay element configurable to impart a predetermined delay to the signal propagating down the delay line. 



2 



EP1 176 439 A1 



Couplers split an optical signal, propagating on the input waveguide, among the tapped delay lines. The delays imparted 
to the signals on each line are chosen so as to effect the desired change to the spectral profile. 
[0010] The claimed embodiments of the invention provide spectral shaping which is polarization-independent, and 
which, in certain embodiments, may be actively configurable or permanently configured. The spectral filtering of the 
invention may, if desired, be provided in a planar device, which is compact, reliable and economically fabricated. It is 
particularly well suited to use as a gain-flattening filter for an optical amplifier, which allows the amplifier to be used in 
a constant gain mode. 

[0011] Additional features and advantages of the invention will be set forth in the detailed description which follows, 
and in part will be readily apparent to those skilled in the art from that description or recognized by practicing the 
invention as described herein, including the detailed description which follows, the claims, as well as the appended 
drawings. 

[0012] It is to be understood that -both the foregoing general description and the following detailed description are 
merely exemplary of the invention, and are intended to provide an overview or framework for understanding the nature 
and character of the invention as it is claimed. The accompanying drawings are included to provide a further under- 
standing of the invention, and are incorporated in and constitute a part of this specification. The drawings illustrate 
various embodiments of the invention, and together with the description serve to explain the principles and operation 
of the invention. Additional features and advantages of the invention will be set forth in the detailed description which 
follows, and in part will be/eadily apparent to those skilled in the art from that description or recognized by practicing 
the invention as described herein, including the detailed description which follows, the claims, as well as the appended 
drawings. 

[0013] It is to be understood that both the foregoing general description and the following detailed description are 
merely exemplary of the invention, and are intended to provide an overview or framework for understanding the nature 
and character of the invention as it is claimed. The accompanying drawings are included to provide a further under- 
standing of the invention, and are incorporated in and constitute a part of this specification. The drawings illustrate 
various embodiments of the invention, and together with the description serve to explain the principles and operation 
of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0014] 

Figure 1 is a schematic diagram of an embodiment of the invention; 
Figure 2 is a schematic diagram of an alternate embodiment of the invention; 
Figure 3 is a schematic diagram of another alternate embodiment of the invention; 
Figure 4 is a schematic diagram of another alternate embodiment of the invention; 
Figure 5 is a graph showing a gain spectrum for an EDFA before gain flattening; and 
Figure 6 is a graph showing a gain spectrum for an EDFA after gain flatteneing. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0015] Reference will now be made in detail to the present preferred embodiments of the invention, examples of 
which are illustrated in the accompanying drawings. Wherever possible, the same reference numbers will be used 
throughout the drawings to refer to the same or like parts. An exemplary embodiment of an apparatus for changing 
the spectral profile of an optical system of the present invention is shown in Figure 1 , and is designated generally 
throughout by reference numeral 10. 

[0016] In accordance with the invention, an embodiment of the invention provides spectral shaping by employing an 
optical tapped delay line realization of the desired narrow-band impulse response. In this embodiment, a splitter tree 
partitions the incoming signal among several taps in the correct proportions for a given impulse response. Each tap 
imparts a nominal delay of T seconds, relative to neighboring taps. The design of the splitter tree may be chosen so 
that the inverse of the delay (1/T) is on the order of the free spectral range of interest, or more preferably, about twice 
the desired free spectral range. Delay perturbations in each tap set the approximate phase for that tap. With the splitter 
ratio determining amplitude, and the delay perturbations determining phase for the signal propagating down each tap, 
the signals output from each tap are superimposed to yield an output signal having the desired spectral shape. 
[0017] As illustrated in Figure 1, an embodiment of the present invention for changing the spectral profile of an optical 
signal includes an array of optical paths 14 which, in the illustrated embodiment, are configured in a splitter tree ar- 
rangement. The optical path array 1 4 can be an array of waveguides in a planar substrate, such as, for example, single 
mode silica waveguides on a silicon substrate, or waveguides on a LiNi0 3 substrate. In the illustrated splitter tree 
arrangement, a first optical path 1 6 directs the signal to a first coupler 1 8, which may be, for example, a Mach-Zehnder 
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interferometer or other suitable variable coupler so that the coupling ratio may be selectively configured as desired. 
The first coupler 18 splits the signal into two portions and couples the two portions into optical paths 21 , 22, each of 
which is optically connected to one of a second and a third coupler 19, 20 respectively, which may be any suitable 
variable 2X 2 coupler such as a Mach-Zehnder interferometer. Each of the second and third couplers 19, 20 split the 

5 power input to it on to a pair of waveguides 24, 26 and 28, 30, respectively. In the illustrated embodiment, these 
waveguides act as taps, with the relative tap weights being effectively set by splitting the input power in appropriate 
ratios over the two stages of variable couplers made up by the first coupler 1 8 (in the first stage) and the second and 
third couplers (in the second stage). Each of the tap waveguides 24, 26, 28 and 30 has a different effective optical 
path length chosen so that each tap imports a nominal delay Tto a signal propagating through it, relative to a neighboring 

10 tap. As illustrated, this may be accomplished by the tap waveguides having successively longer (or shorter) actual 
lengths. For purposes of illustration, a two-stage splitter tree is illustrated in Figure 1. Embodiments with different 
numbers of stages and of tap waveguides are within the scope of the invention. 

[0018] A delay element 38 is associated with each of the tap waveguides 24, 26, 28 and 30 for imparting delay 
perturbations to the signal traveling through that tap. The delay elements 38 may be, for example, portions of the 

15 associated waveguide which have an altered index of refraction. Altering the index of refraction in a waveguide changes 
the effective optical path length through that waveguide, thereby effecting a delay in propagation time, and in phase, 
through that waveguide. The alteration in refractive index can be a permanent alteration, or alternatively, a temporary 
or reversible alteration. In this alternative, the refractive index preferably could be selectively modified as desired to 
achieve the desired delay in propagation. 

20 [0019] Some examples of how this modification of refractive index can be achieved include, for waveguides in a 
silica substrate, exposing a portion of the waveguide to UV radiation. This will result in a permanent modification of 
the refractive index. Temporary or reversible modifications to the refractive index can be provided for waveguides in a 
silica substrate by providing selectively actuated heating elements adjacent the waveguide portion for which the re- 
fractive index changes are desired. For an embodiment of the invention using waveguides in an electro-optic substrate 

25 such as UNi0 3 , electrodes adjacent a waveguide may be used to provide a voltage across the waveguide, which will 
change the refractive index in a non-permanent fashion according to the well-known electro-optic effect. 
[0020] In the illustrated embodiment, the tap waveguides are reflected back on themselves by, for example, a refledtor 
34 at the output end of the optical paths which make up tap waveguides 24, 26, 28 and 30, reflect the portions of the 
signal propagating down these optical paths back through the splitter tree 14 in the opposite direction. This ensures 

30 that the couplers, which separate the input signals into the-successive stages of the splitter tree and the coupling of 
the delayed copies of the input signals, are accomplished in the same ratios. The reflector 34 can be a polished gold- 
plated surface or any other suitable reflective optical element or elements. The ports 40, 42 and 44 to the first, second 
and third couplers 18, 19 and 20, which do receive input power may be used to probe the couplers, for system moni- 
toring. 

35 [0021] Figure 2 shows an embodiment of the invention similar to that of Figure 1 , in which the first, second and third 
couplers 18, 19, 20 are illustrated as Mach-Zehnder interferometers 41 , 43, 45, respectively. In an embodiment where 
the splitter tree 14 is fabricated with silica waveguides on a silicon substrate, the splitting ratios through these Mach- 
Zehnder elements 41, 43, 45 may be modified as desired by actuation of heating elements 46, 47, 48 adjacent the 
waveguide making up one arm of the Mach-Zehnder interferometer. In the embodiment of Figure 2, a circulator 50 

40 between the input waveguide 52 and the splitter tree 14 ensures that the reflected signal propagates to the output 
optical path 54. Those skilled in the art will appreciate that alternative means may be used to direct optical signals from 
a communication system into and out of the device 10, employing other optical components such as isolators, for 
example. 

[0022] The embodiments of Figures 1 and 2 may be modified by replacing the reflector element 34 with a mirror 

45 image of the array. An example of this is illustrated in Figure 3. 

[0023] In another aspect of the invention, illustrated in Figure 4, an optical communication system includes a trans- 
mitter 53,- a receiver 54, an optical amplifier 56 which may be an erbium-doped fiber amplifier (ED FA), and a gain- 
flattening" filter 58. The gain-flattening filter can be a spectral filter 1 0 as illustrated in Figures 1 , 2 or 3. A waveguide 
60 propagates a signal from the amplifier 56 to the gain flattening filter 58. In this aspect, the splitting ratios and delays 

so for the embodiment illustrated In Figure 1 , 2 or 3 would preferably be chosen to provide a gain flattening effect to the 
output of the optical amplifier 58. 

[0024] For example, suppose it is desired to realize a gain-flattening filter for an erbium-doped fiber amplifier (ED FA), 
operating over the C-band (roughly 1530-1560 nm), with 8 taps. In this particular example, the parameters will be 
chosen so that the device optimally flattens the gain of the amplifier at a population inversion level of 65%, the expected 
55 operating condition for a particular amplifier design. Figure 5 shows the spectral gain profile for the amplifier at that 
inversion level. It is recognized that the spectral measurement is noisy, and the optimization should be robust enough 
to accommodate this limitation. 

[0025] Calculation of the appropriate complex tap weights may be done using MathCAD software, and an annotated 
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MathCAD script is given in Appendix - Part 1 to illustrate an example of appropriate steps taken in designing a device 
to realize such a function. Spectral scans at several inversion levels are recorded in the file, C-Band.txt, from which 
the design wavelengths (indicated by the vector k - Design) and gains (indicated by the matrices: y in dB units, A in 
absolute units) are extracted. 

5 [0026] The Nyquist (sampling) theorem ordinarily dictates that the sampling interval T be chosen to correspond to 
the inverse of twice the frequency range of interest. In the present design example, this criterion would lead to the 
specification T = 110 fsec. However, this should not be interpreted as a rigid constraint, and in fact a minimum least 
square error solution is actually obtained from 1 85 fsec. In fact, if the Nyquist theorem is followed blindly, there is also 
the very real possibility that a quite lossy system will result, since the range over which the desired response is specified 

*o is only half the free spectral range (in fact, it is this condition which insures that there is absolutely no aliasing), and 
over the other half may be where the transmission maximum occurs. This possibility follows from the constraint that a 
given 'tap weight cannot exceed unity. Thus, some searching for an optimal T must occur in general. For a given T 
however, optimization proceeds as follows. 

[0027] The first step is to quickly arrive at a linear least mean-square (LMS) solution. This leads to a solution for the 

15 complex tap coefficients which is obtained by representing the system transfer function as a weighted (by the tap 
coefficients) superposition of sinusoids in the frequency domain. Under such an optimization paradigm, the transfer 
function to be approximated must be specified both in amplitude and phase. Since the design task simply specifies 
the amplitude of the desired spectral response, the phase can be any arbitrary function. The phase functions <D in 
Appendix — Part 1 were chosen based on a minimum-phase criterion. In the same way that the real and imaginary 

20 parts of the refractive index are related in dielectric materials, the principle of causality leads to this relationship between 
the amplitude and phase of a linear, causal system. Clearly, this is not the only phase relationship that can exist for a 
given amplitude of response, and in fact a non-linear optimization leads to a different function, as will be explained 
shortly. However, it is necessary to initiate the optimization at a reasonable starting point, and in this example, this 
function provides that initial guess. 

25 [0028] Thus, with amplitude and phase determined, the target transfer function H is completely specified. The graphs 
of ^design k versus IH k 3 I and 4 K 3 illustrate the amplitude and phase for this function, chosen for N=3, where the 
population inversion is 65 %. In the present example, the number of taps, indicated by the constant, Taps, is 8. The 
complex cumulants are then computed for the LMS estimation of the complex tap weights, represented by the vector 
h, which is then normalized and presented graphically in the MathCAD script. From this, two frequency responses are 

30 derived: H2 and H3. H3 represents the response of an ideal tapped delay line system, for which h is not frequency 
dependent, while H2 represents the response for a system in which the phase is allowed to vary over the frequency 
range in a way consistent with approximating a desired phase by a constant delay over a narrow range of frequencies. 
The graphs of A._Data m verses I H 2 m l and I H 3 J show the difference between the two functions. There is little difference 
between the functions, validating the assumption that for this narrow-band system, a constant delay is a good approx- 

35 imation to a constant phase. Over the frequency range of interest in the example, the RMS ripple in the compensated 
output is only 0.2 dB. 

[0029] As mentioned above, the linear prediction may be used as a starting point for optimization for a given design. 
Once h is found in this way, a non-linear optimal solution may be found by entering the elements of h into the C program 
shown in Appendix - Part 2, which randomly perturbs the real and imaginary components of h in such a way as to 
40 minimize the RMS error in the gain-compensated ripple, after an appropriate number of trials and with simulated an- 
nealing. 

[0030] Following this non-linear optimization, the optimal impulse response coefficients (i.e. the elements of h), given 
in Appendix-Part 3, show that the RMS ripple has been reduced to a mere 0.047 dB. As can be seen in Figure 6, this 
ripple is most likely due to random error in the data. 
45 [0031] It will be apparent to those skilled in the art that various modifications and variations can be made to the 
present invention without departing from the spirit and scope of the invention. Thus, it is intended that the present 
invention covers the modifications and variations of this invention provided they come within the scope of the appended 
claims and their equivalents. 

50 
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Appendix - Part 1 

Data which, in this example, consists of gains for 71 wavelengths at 11 inversion levels, is read in: 
Spectra 5"C-Band.txT 
V : s READPRN( Spectra) 

M:=71 N = 11 

m := 0.. M - 1 
n:=0 r N-1 
Data wavelengths 

X Data m := V_ * «. 0 

Inversion levels 

n2 n :=V 0,n^15 

Design wavelengths. In this example, only the useful portion of the gain curves are of interest, 
k = 0.. 34 

X_Design k := *<_Data k 
Radian frequencies: 

2-K-3-10 8 

m X Data m .10-° 
— m 

Cains at k'th wavelength, nth inversion level, in dB 

r k,a ,sV k+1,nt 15 

Corresponding sampling rate, in seconds. 

T~:= 1.85.10 - 13 
Absolute spectral response. 

Tk.n 

\.n"« 10 
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Mean radian frequency for n'th inversion level. 



10 



k 



Create a phase profile, based on a minimum-phase criterion. Assume spectrum to repeat after a free 
spectral range 



15 



20 



25 



30 



40 



5 := 0.. 34 p =-5-5 



k.n 

P 



(p-35+.k-s)»0,0, 



p-35 + k-s. 



♦k.n := mod («k.n'«- 2 ) 



Desired transfer function for gain compensation at n'th inversion level (choice of numerator sets desire< 
gain). 

■K.n 
0.06 



H kf3 [0.04h \ 



0.02 



"i r 




0.5 



♦ k.3 o 



J I I L 



1520 1530 1540 15SO 15C 



Number of waveguide taps (in this example, 8 taps). 
Taps:= 8 



-0.5 



"i 1 1 r 




J I I L 



1520 1530 1540 1550 1560 1 
l_Oeslgnfc 
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10 



15 



20 



25 



v— 

30 1 Norm 



35 



40 



50 



55 



Choice for an inversion level at which to operate, 
N := 3 
n2 N = 0.65 

Frequency response vector and matrix for least mean square (LMS) error computation, 
t := 0 Taps - 1 u := 0 .. Taps - 1 

k 

k 

Impulse response vector which gives the LMS error, 
h := IVT 1 .v 



Normalize impulse response. 

Norm:=£|h t | 
t 



Amp!itude t := | h t | 



Amplitude^]) 1 »1 

t 



JDOCID: <EP 
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Amplitude and phase of the impulse response which gives the minimum RMS error for the 
chosen number of taps, at the chosen inversion level: 



10 



15 



20 



25 



h t | 0.5 



b 5 
0.664294 + 4.84072 -10~ 3 i 

-0.041764+ O.04456i 

.0.1 07477- 0.061 672i 
0.0624+ 0.076113i 

-0.018686+ O.G23785i 

-4.532982-10" 3 - 4.09261 5- 1 0 r3 I 
-4.763722-1 0~ 3 + 0.0125961 
-7.402329* 10* 3 - 4.12571 1 *10"* 3 i 



ar 0 (h t ) 




30 



Actual frequency response of Tapped Delay Line.(TDL) system. 



H v s £l h tl- e *p 



•«p(-]-» m -tT) 



40 



Frequency response of ideal LMS TDL system, 
t 



45 



50 



I H2 ni I 0.5 



i r 




\ I 

I H3 m I 0.8 Y 'J 



1520 1540 1560 1580 
*-Data m 



1 520 1 540 1 560 1 580 160 
X_Data m 
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Gain compensation: Amplifier + TDL filter: 



Comp k :=10.|og[|A ktN .(]H2 k p 2 |] 




Normalize Impulse response. 



Norm:= £(] h tD 
t 





■ A. 

Norm 




arg(M 


0.664311 


7.286887-10° 


0. 061072 


2.323812 


0.123914 


-0.520946 


0.092406 


0.967871 


0.030248 


2.236701 


6.107161 -10* 


-2.407204 


0.01346 


1.932358 


8.47443 -10* 3 


-2.633122 - 



ji - 24.340649 



max(Comp) - mln(Comp) ■ 0.675054rnax Ripple indB 



2(Comp k - M ) 5 
k 



0.207617 RMS Ripple 



10 



EP1 176 439 A1 



Appendix - Part 2 

/• AGPf.l.e ♦/ 

/• This 'program performs non-linear optimization of the complex tap weights 
for the active gain flattening filter, starting from the LMS predicted 
tap weights from the MathCad simulations. 
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tine lude <stdio.h> 
# include <stdlib.h> 
I include <math.h> 



/• Constants */ 



int const Taps - 8; 
double const T - 1.85E-13; 
15 int const Trials • 10 000; 

double const Step - 0.001; 



/* Sampling interval. 

/* Number of trials in search for optimal tap coefficients. 
/* Step size for coefficient perturbation. 



20 



25 



30 



/* Variables */ 

int ra,n,k, Laatjnj 
float Dummy j 
int Num; 

int Lambda (100) i 
double Onega [100] \ 
double argtioo] ; 
double Gain [100] ; 
double Target [100] ; 
double Hurdle; 
double h_Ref ; 
FILE 'Spectrum; 
FILE •Outtile; 
double h(l6] - ( 



0 .968829, 
-0. 060909, 
0,156748, 
0.076421, 
-0.027252, 
-0. 00661106, 
-0.0069476, 
-0.010796, 



/♦ Number of frequencies considered. 

/« Wavelengths in nm. Contained in first column of the file Spec 

/* Corresponding omega's. 

/* Product of omega and sampling interval. 

/* Successively lowest XLEP values. 

/• Dummy variable for perturbing Re and Xm parts of impulse reap 

/• File contains spectrum to be flattened. 

/* File contains optimized impulse response; 

0.00705968, 

0.064988, 

-0.089945, 

0.111006, 

0.03469, 

-0.00596BB13, 

0.016371, 

-0.006017082); 



35 



40 



45 



50 



/* Functions */ 

double SUnt ComputedAt); 

double ILEFt) ; 

main(argc, argv) 
int argc; 
char «»argv; 

{ 



/* Computes the spectrum for the as* then defined 

complex impulse response, at the radian frequency 
indexed by the integer ComputedAt. 

/♦ Resultant ILEP. 



Spectrum • fopen(»a:\Gain. txt - , # r*) / 
Outfile - fopen fa: \h.txt", •*■) ; 



/* Read in gain data. First column is wavelength, second is gain, in dfl. 

The wavelengths should only cover the range we're interested in flattening. 
After reaching EOF, the variable Num will represent the number of wavelengths 
represented in this range. 



Num • 0; 

while (fscanf (Spectrum, "%d %f \n\ fcLambda (Num] , tDumray) !■ 

Gain (Num] - pow(l0, Dummy/10} ; 

Omega (Num) » 1. 884956E18/ (float) Lambda [Num] ; 



EOF) 
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arg (Numl - Omega [ Mum] * T ; 
Num>*; 

/* compute inverse function for 25 da overall gain •/ 
for(m-0;m<Num;m+t) Target(m) - 316/Gain(ra]j 

Hurdle • ILEP()r /* Initiates search. */ 

I*ast_m • 0;m • 0; 
for (n-0;n<Triala;n*«.) 

' while (m— Last_m) m - (int) floor (2 . 0*Tapa*rand {) / (float) RAND_J4AX) ; 
h Ref - Mini ;~ 



him) V- Step; 
while(ILEP()<Hurdle) { 



h[m) ■ hjtlef;. 
„h[n) Itep; 
while (ILEP() <Hurdle} { 



him) • h_Ref; 
Last m - m; 



Hurdle - ILEPO ; 
h Ref - htm) ; 
hTm) Step; 



Hurdle = ILEPO ; 
h_Ref ■ h(m] ; 
hTm) Step; 



} 

printfC « \n», Hurdle) ; 

for(m-0im<Taps/m**) printfC %f %f \n« ,h|2»m) , h(2*ro*l) ) ; 
for(m-0;m<Taps;m^) fprintf (Outf ile, ■ %f *f \n«,h(2*m) ,h[2«m+l) ) ; 

f close (Spectrum) ; 
f close (Outf ile) ; 

} 

double Stint At) 
{ 

double Dummy _Re, Duminy_Ino; 

Dumray_Re • 0? 

Dummy_Im ■ 0; 

for (k-0;k<Taps;k+-f) 

Dummy Re h(2*k) ♦cos(arg[Atl *k) ♦ h (2*k*l) *sin(arg (At) *k) ; 
Dumray'lm h [2«k*l) *cos (arg (At) *k> - h(2*k) *sin(arg(At] «k) ; 

return (Dummy_Re*Dummy_Re Duireny_Im*Duramy_Im) ; 

double ILEF (> 
{ 

double min, max, Dummy; 
int m; 

min - 10*logl0(S(0)«Gain[0]); 

max ■ min; 

f or ( m- 1 ; m<Num ; m+ + ) 

Dummy . 10*loglO (S (ra) *Qain (m) ) ; 
if ( Dummy > max) max « Dummy; 
if (Dummy <min) min » Dummy; 



return (max -rain) ; 
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Appendix- Part 3 

Amplitude and phase of the impulse response which gives the minimum RMS error for the 
chosen number of taps, at the chosen inversion level: 
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25 



M 0.5 



0 5 

0.654353+ 9.581827-10" 3 i 
-0.03699 + 0.038666i 

0.115147 -0.071 02i 

0.054962 + 0.088909i 
-0.020698+ 0.019021 

-2.750918»10" 3 + 1. 036431 •10~ 3 i 

"7.045228-10"* 3 + 0.01 1 345t 



-7.764422-10 



•Iff 3 - 



1. 167252*1 0" 3 ! 
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Actual frequency response of TDL system. 
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Frequency response of ideal LMS TDL system. 
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Gain compensation: 
Amplifier + TDL filter: 
Comp k :=lO.|og[|A k(N .(]H2 k D 2 |] 



E C ° mp k 
k 



Normalize impulse response. 

Norm:=£(]h t D 
t 



Ht Norm 



Ih, I 


arg( h, ) 


0.654423 


0.014642 


0.05351 


2.334039 


0.135287 


-0.552664 


0.104525 


1.017118 


0.02811 


2.398415 


2.939684-10"" 


2.781281 


0.01335 


2.126527 


7.85167M0" 0 


-2992377 



u,= 24.169417 
max(Comp) - min(Comp) ■ 0.1 5693&nax Ripple in dB 



Z( CORiP K^) 2 

k 



E 1 

k 



: 0.047031 RMS Ripple 



Claims 

1. An apparatus for changing the spectral profile of an optical signal comprising: 

a plurality of optical paths, each of said optical paths having associated therewith a delay element which 
imparts a predetermined delay to a signal propagating through the optical path; 

a plurality of couplers configured to split the optical signal among said plurality of optical paths according to 
predetermined splitting ratios; 

wherein said predetermined delays and predetermined splitting ratios are chosen so as to effect the desired 
change to said spectral profile. 

2. The apparatus of claim 1 wherein said predetermined delays and said predetermined splitting ratios are variable. 

3. The apparatus of claim 1 wherein said optical paths comprise waveguides and said delay elements include portions 
of said optical paths with an altered index of refraction. 

4. The apparatus of claim 3 wherein said altered Index of refraction is permanently altered. 

5. The apparatus of claim 4 wherein said waveguides comprise silica waveguides, and said portions of said optical 
paths with an altered index of refraction comprise portions of said silica waveguides which have been exposed to 
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UV radiation. 

6. The apparatus of claim 2, wherein said waveguides comprise waveguides in electro-optic material and said delay 
elements include electrodes adjacent said waveguides. 

5 

7. The apparatus of claim 3, wherein said waveguides comprise silica waveguides and said delay elements include 
heating elements adjacent said waveguides. 

8. The apparatus of claim 3 wherein said index of refraction is.selectively reversibly Talterable. 

10 

9. The apparatus of claim 3, wherein said optical paths comprise waveguides forming a splitter tree. 

10. The apparatus of claim 8 wherein said optical paths comprise waveguides forming a splitter tree. 

is 11. The apparatus of claim 10 wherein said splitter tree has an input end and an output end, and further comprising 
a reflector at said output end. 

12. The apparatus of claim 9 wherein said splitter tree has an input end and an output end, and further comprising a 
reflector at said output end. 

20 

13. The apparatus of claim 1 further comprising an optical amplifier and an optical path between said optical amplifier 
and said coupler. 

14. A gain flattening filter for use with an optical amplifier comprising: 

25 

an input waveguide; 

a plurality of tapped delay lines comprising waveguides and delay elements associated with said waveguides, 
a plurality of variable couplers configured to split an optical signal propagating on said input waveguide among 
said plurality of tapped delay lines according to predetermined splitting ratios; 

30 . 

wherein each of said delay elements is configurable to impart a predetermined delay to the signal propagating 
down the associated waveguide, wherein said predetermined delays are chosen so as to effect a flattening of the 
spectral profile of said optical signal. 

35 15. The gain-flattening filter of claim 14 wherein said tapped delay lines comprise waveguides formed in a planar 
substrate. 

16. The gain-flattening filter of claim 14 wherein said delay elements comprise selectively actuated resistive heating 
elements near said associated waveguides. 

40 

17. The gain-flattening filter of claim 1 4 wherein said delay elements comprise regions of said associated waveguides 
having a propagation coefficient which differs from that of the remainder of said associated waveguide. 

18. The gain-flattening filter of claim 16 wherein said regions of said associated waveguides comprise regions which 
45 have been exposed to UV radiation. 

19. An optical communication system comprising: 

an optical amplifier; 
so a spectral filter comprising an input waveguide; 

a plurality of tapped delay lines, each of said tapped delay lines including a delay element on said line, 

a plurality of couplers configured to split an optical signal propagating on said input waveguide among said 

plurality of tapped delay lines according to predetermined splitting ratios; 

55 wherein said delay elements are configurable to impart a predetermined delay to the signal propagating down said 

line, wherein said predetermined delays are chosen so as to effect the desired change to said spectral profile; and 
an optical path between said amplifier and said filter. 



15 



EP1 176 439 A1 



20. A method for modifying the spectral distribution of an optical signal comprising the steps of: 

splitting the signal among a plurality of optical paths, each having an output end, whereby a portion of the 
signal propagates on each of the optical paths; 

changing the relative phases of the portions of the signal propagating through the optical paths; 
modifying the amplitudes of the portions of the signal propagating through the optical paths; 
the relative phase changes and the amplitude modifications having magnitudes which result in a desired spec- 
tral distribution when the portions of the signal at the output ends of the optical paths are combined; and 
combining the portions of the signal at the output ends of the optical paths. v - 

21. T,he method of claim 20 wherein the step of changing the relative phases comprises delaying the propagation of 
the portion of the signal in at least one of the optical paths by an amount different from the others of said optical 
paths. 

22. The method of claim 21 wherein the amplitude modifying step comprises applying weighting coefficients to the 
portions of the signal propagating through the optical paths. 

23. The method of claim.22 wherein the delaying step further comprises modifying the refractive indices of each of 
said plurality of the optical paths. 

24. The method of claim 23 further comprising the steps of: 



determining the magnitude of each of said delays necessary to result in the desired spectral distribution; and 
modifying the refractive index of each of the optical paths by an amount which results in the determined delay 
magnitude. 



25. A method for modifying the spectral distribution of an optical signal comprising the steps of: 

splitting the signal among a plurality of optical paths having output ends according to predetermined splitting 
ratios, 

changing the relative phases of the portions of the signal propagating through different ones of the optical paths; 
the relative phase changes and the splitting ratios having magnitudes which result in a desired spectral dis- 
tribution when the portions of the signal at the output ends of the optical paths are combined; and 
combining the portions of the signal at the output ends of the optical paths. 
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FIG. 1 
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i SDOCID: <EP 1 176439 A 1_L> 
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FIG. 4 
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